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A/'ms. To use a 3D numerical MHD experiment representing magnetic flux emerging into an open field region as a background field for 
tracing charged particles. The interaction between the two flux systems generates a localised current sheet where MHD reconnection 
takes place. We investigate how efficiently the reconnection region accelerates charged particles and what kind of energy distribution 
they acquire. 

Methods. The particle tracing is done numerically using the Guiding Center Approximation on individual data sets from the numerical 
MHD experiment. 

Results. We derive particle and implied photon distribution functions having power law forms, and look at the impact patterns of 
particles hitting the photosphere. We find that particles reach energies far in excess of those seen in observations of solar flares. 
However the structure of the impact region in the photosphere gives a good representation of the topological structure of the magnetic 
field. 
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1. Introduction 

To investigate the dynamical evolution of the large scale phe- 
nomena seen in the solar atmosphere, one traditionally adopts 
the magneto-hydrodynamic (MHD) representation. The MHD 
approach is a macroscopic approximation describing the time 
dependent evolution of a plasma. The assumption when adopt- 
ing MHD is that the time dependent plasma evolution is well 
described by macroscopic parameters, such as temperature, den- 
sity, bulk velocity and magnetic field. This requirement is ful- 
filled when the plasma is in a local thermodynamic equilibrium 
and the characteristic MHD length scale is much longer than the 
mean free path of the plasma particles. For most of the time the 
plasma in the lower solar atmosphere obeys this requirement, but 
there are special situations for which MHD is not a good repre- 
sentation. One such situation is magnetic reconnection. Despite 
this fact, magnetic reconnection has for many d ecades been in- 
vestig ated using the MHD approximation (See Priest & Forbes 
2000, for a review.). In MHD the non-ideal evolution imposed 
by viscosity and resistivity in the plasma is described through 
a parameterization that approximates the real particle interac- 
tion on length scales well below the typical resolution of MHD. 
The choice of the form and value of the magnetic resistivity, 77, 
naturally has implications on how the reconnection process pro- 
gresses. 

From the vast literature on 2D MHD reconnection, a com- 
mon feature is the presence of a localised diffusion region, from 
which standing slow mode sho cks extend alon g the separator 
lines - the Petschek like picture (Petschek 1964). When redoing 
this type of experiments with particle codes, one finds a differ- 
ent picture. The existence of the MHD slow shocks vanishes as 
the evolution of protons and electrons decouples on sufficiently 
short length scales removing the plasma's ability to make a sin- 



gle coherent shock structure dBirn & Pri est 2007, See the article 
by Drake and Shay). Particle simulations allow for distribution 
functions to reach significant deviations from Maxwell distribu- 
tions, typically providing a power law tail of high energy parti- 
cles accelerated away from the diffusion region. One could there- 
fore fear that the global evolution of the field will be qualitatively 
different depending on which method one chooses to use for the 
investigation. Recent comparisons between different approxima- 
tions, going from MHD to particle in cell (PIC), in a controlled 
reconnection experiment have shown that despite variations in 
the reconnection speed between the different approaches, the 
same amount of flux reconnects a nd the same gene ral structure 
of the final magnetic field results (Birn et al1 l2005l) . Numerical 
MHD experiments therefore do capture the large scale dynamics 
of the evolution, and will at worst provide a somewhat wrong 
time scale and energetics for the dynamical evolution. 

In the Solar atmosphere magnetic reconnection is deemed 
responsible for various explosive events, and the conversion of 
magnetic energy into other forms is therefore an important pro- 
cess one needs to understand to be able to explain observations. 
In the MHD approximation magnetic energy is converted to bulk 
motion and Joule heating of the plasma. This is a limitation, 
that already became clear from observations of X-rays and y- 
rays in the late 1970's. Already then it was found that a signif- 
icant amount of the energy release associated w ith large flares 
is converted into ene rgetic particles (10-50%) dLin & Hu dson 
119761: iLin et al.ll2003l) . The effect of the accelerated particles is 
seen in different characteristic features. In the general picture 
one assumes that reconnection takes place in the corona, and 
that the strong electric field parallel to the magnetic field ac- 
celerates electrons and protons to velocities where their cross 
section becomes so small that they interact very little with the 
coronal plasma. Only when they reach the dense plasma close 
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to the photosphere are they stopped very efficiently while, caus- 
ing them to emit bremsstrahlung radiation. In loop flares, this 
gives rise to H a ribbons in the photosphere that move apart with 
time. The combined information of the motion and the mag- 
netic field distribution in the photosphere provides indirect infor- 
mation about the reconnection speed in the acceleration region 
dHovng et alJll981t l lsobe et alJl2002t Sing et alJl2007t fLiu et all 
20Q8ar iLiu & Wangll2009l) . In the reconnection model, one also 
expects the particles to be accelerated along the open field lines 
reaching into the upper corona. H ard X-ray observat ions have 
been made at the loop top of flares (Masuda et al. 1994]), but it is 
not until recently that RHESSI has been able to detect a signature 
of independent particle distribution s on either side of an assumed 
current sheet dKrucker et alj|2008h . This investigation strongly 
supports the picture where the acceleration region is located in 
the corona, and that particles are accelerated in both directions 
away from the local reconnection region. Their analysis also in- 
dicate that the very energetic particles can be maintained in the 
local vicinity of the diffusion region long enough to be able to 
interact with the local plasma before continuing away from this 
region. Understanding the particle acceleration and interaction 
process is therefore important for understanding the distribution 
of the released magnetic energy in the flare process. 

The lack of information regarding particle acceleration in nu- 
merical MHD experiments has led us to start investigating its 
effect. We adopt a simplified approach, where test particles are 
traced in snapshots from a numerical MHD experiment. The aim 
is to investigate the efficiency of the acceleration through form- 
ing distribution functions, looking at the patterns they form as 
they exit the allowed tracing domain and especially calculate ra- 
diation spectra from the bremsstrahlung that arises when the fast 
electrons are stopped in the top layer of the dense photosphere. 

The layout of the paper is as follows. In the following section 
we review the relevant information about the numerical MHD 
experiment used as background field for test particle tracing. 
Section [3] describes the numerical approach used for tracing the 
particles, while Section [4] gives the setup of the tracing experi- 
ment. In Section|3J the results of the data analysis are presented. 
Section [6] discusses the drawbacks with the approach and what 
we need to do to improve the results, while Section [7] draws up 
the conclusions of the paper and suggests further work to be 
done. 

2. Numerical MHD experiment 

The dynamical evolution of an emerging flux tube entering into 
an open coronal fi eld was investigated u s ing a numerical 3D 
MHD approach by [Moreno-Insertis et alj d2008l) . The numeri- 
cal experiment consists of a stratified atmosphere covering the 
region from 4 Mm below the photosphere to 29 Mm above it. 
Imbedded in the convection zone is a twisted flux tube that over 
a fraction of its length is made buoyant. Additionally a slanted 
magnetic field permeates the numerical domain, representing an 
open magnetic field configuration. The part of the flux tube that 
rises from the convection zone pushes away the slanted magnetic 
field to make space for itself. While the tube is still in the convec- 
tion zone, the motions are slow and not very compressible, and 
very little interaction between the two flux systems takes place. 
As the rising tube reaches the photosphere the magnetic field 
strength starts to build up. Shortly after, the vertical magnetic 
pressure gradient becomes sufficient to initiate the emergence of 
the flux tube into the upper atmosphere. The emergence speed in 
this phase is much higher than the slow rise speed observed in 
the convection zones. The initial plasma ft decreases with height 



through the transition region and into the corona making it eas- 
ier for the two magnetic flux systems to come into close contact. 
The change in f3 is important for the dynamical evolution of the 
boundary layer between the two flux systems, especially where 
they experience large differences in field line orientation. Such 
regions are prone to generating strong current concentrations as 
they are compressed by the emergence process. As part of the 
flux tube emerges into the corona a current dome separating the 
two flux systems is formed. The continued emergence concen- 
trates the current in the dome into a long thin sheet connecting 
along the length of the emerged flux, laying skew to one side of 
the summit line of the emerged flux. Eventually the current den- 
sity becomes large enough for non-ideal MHD processes to start 
changing the field line connectivity. The reconnection between 
the magnetic field from the rising loop and the ambient open 
field creates a new low lying loop system next to the emerg- 
ing one and new field lines that connect the deeper parts of the 
twisted loop with the open flux. The effect of reconnection is 
therefore to eat away the emerging flux and build up a secondary 
loop system. One way to see this is to investigate how the sepa- 
ratrix surface between the two flux systems changes with time. 
This is investigated by Galsgaard & Moreno-Insertis (2010, in 
prep.) and shows how the interaction between the magnetic field 
and the photosphere changes structurally with time, clearly in- 
dicating the continued diminishing of the area occupied by the 
emerging loop system and the concurrent growth of the new loop 
system. Fig. Q] shows the connectivity domains at three differ- 
ent times during the emergence process. The background shad- 
ing represents the photospheric magnetogram indicating the sign 
and magnitude of the magnetic fields normal component. This 
shows a bipolar structure imbedded in a monopolar background 
field. From the magnetogram alone one could be led to conclude 
that a new bipolar field has emerged into the photosphere with 
no complicated connectivity pattern. This is not the case: The 
lines on the map show the intersection of the separatrix surfaces 
with the photosphere. The main features are the outer line defin- 
ing the perimeter inclosing flux that connects both ends to the 
photosphere. A line dividing the perimeter represents the divi- 
sion between the two loop systems, where the top one is the sec- 
ondary loop created by the magnetic reconnection process. The 
jagged form of the outer line is created by convective motions 
that overshoot above the plane in which this surface is shown. 

As a consequence of the changes in magnetic connectivity, 
two reconnection jets are accelerated out of the current sheet 
edges. The most prominent of the jets escapes upwards, giving 
rise to high velocity plasma motions that follow the open field 
lines away from the solar surface. 

The evolution of the numerical MHD experiment resembles 
the anemone jets which Hinode o bservations have fo und to be 
very common in open flux regions (Shibata et al. 2003). 

The MHD approximation assumes the plasma to be in a sta- 
tistical equilibrium, defined by frequent collisions between par- 
ticles, where physical variables like temperature and bulk ve- 
locity are represent for the average particle in a small volume. 
Therefore, numerical MHD experiments are not able to provide 
indications as to how individual particles may react to strong 
electric fields parallel to the magnetic field. In a real plasma, 
this electric field component is responsible for accelerating a 
small fraction of particles up to relativistic velocities. As only a 
small fraction of the particles are accelerated, they give rise to a 
high energy power law tail attached to the Maxwellian velocity 
distribution. Power law distributions created in eruptive events 
like solar flares are well documented by RHE SSI observations 
dLiu et al . 2008bi lLinl2008tlLiu & Wand2009h . In fact, integrat- 
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Fig. 1. The connectivity pattern of the magnetic field in the emer- 
gence process at three different times. The colour indicates the 
field strength and sign of the normal component of the magnetic 
field with blue representing negative flux. Each frame is scaled 
to its own dynamical range. The lines indicate the intersection of 
the separatrix surfaces with the photosphere. 



ing the ene rgy content of this tail indicates that a larg e fraction, 
up to 50% (iLin & Rhessi Teaml2002t ILin et al.ll2003h . of the en- 
ergy released in the flare process may go into accelerating par- 
ticles. If we therefore want to understand in detail the effect of 
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the flare processes it is important to take this additional physics 
into account. This is not an easy task, as it requires models to re- 
solve particle interactions, while still maintaining the large scale 
structures responsible for the dynamical driving of the magnetic 
system. Resolving both the particle scale and the large structural 
scales requires unrealistically large numerical resolution. A sim- 
pler approach is to assume: 

- The integration time for the particle experiment is much 
shorter than the time between MHD snapshots. 

- The MHD snapshot is static in time. 

- There is no back reaction on the MHD field from the particle 
motions. 

The first assumption is fulfilled in the following analysis, while 
the two later assumptions are unrealistic. However, using these 
makes it possible to obtain results regarding (i) the particle dis- 
tributions after acceleration, (ii) the locations where the acceler- 
ation takes place and (iii) where the accelerated particles end up 
as they are rapidly accelerated away from their initial positions. 
The details of the applied approximation for this approach are 
discussed in the following section. 

3. Test particle tracing approach 

The relativistic equation of motion for a charged particle in the 
presence of a Lorentz force is described by the Lorentz equation 
(LE): 

j=?(E + vxB), (1) 
at 

where p = myy is the relativistic momentum, q the particle 
charge, E and B the electric and magnetic field at the particle 
position, v is the particle velocity, m is the rest mass of the par- 
ticle and y is the Lorentz factor. This equation describes the or- 
biting motion of a charged particle around any magnetic field. 
Before going further into the discussion of the applications of 
this equation it is worth looking into the characteristic radius of 
the orbiting motion of a charged particle in a magnetic field. The 
Larmor radius of this circular orbit is 

Rc = ^, (2) 
Bq 

where v is the speed. For the emergence experiment the charac- 
teristic coronal values are B = 10~ 2 T and the coronal temper- 
ature 1.2 MK. Assuming a Maxwell-Bolzman distribution for 
the particles gives typical electron velocity of roughly lQ 7 m/s 
and a y — 1 . Using these values in the equation above, we ob- 
tain a Larmor radius of order 10 3 m. This is a typical value, but 
more extreme circumstances, may change it by up to 2 orders 
of magnitude. The relevant quantity to compare to this value 
is the resolution of the emergence experiment. In this case the 
box size is 34Mm and the grid resolution on the order of 320 
grid points. This gives a typical grid resolution in the numeri- 
cal domain of order O.lMm. Comparing numbers it is seen that 
the orbiting motions typically take place on 10~ 7 of a grid cell. 
Obtaining the magnetic field with this spatial resolution, from a 
single precision dataset, implies that the magnetic field seen by 
a particle in a given orbit, will be constant across the orbit, and 
its orbit will only change due to changes in the magnetic field 
direction taking place over large distances measured in Larmor 
radii. Using Eq. (Q3 to track particle motions in the numerical 
MHD experiment will therefore require a huge amount of com- 
puting time, while providing information that can be obtained 
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using a simpler approximation. Instead of using the Lorentz 
equation to determine the part icle paths, we u se the Guiding 
Center Approximation (GCA) dNorthroplll964l) . The GCA in- 
tegrates the B-parallel motion and the non-parallel of charged 
particles, while ignoring the orbital motion, but keeping track of 
the magnetic moment. Ignoring the orbital motion implies that 
the time step for this approach can be up to several orders of 
magnitude larger than for the Lorentz equations, which can save 
a huge amount of computing time and provide the same particle 
trajectories. The relativistic GCA equations of motion are: 



du\\ 
dt 
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where the y refers to the component along the magnetic field, 
R is the center coordinate of the particle orbit and j_ refers to 
the component perpendicular to the magnetic field vector - here 
only the magnitude is interesting and not it's angle. These equa- 
tions represent the momentum equation, the motion of the guid- 
ing center and the adiabatic invariance of the magnetic moment. 

We solve these GCA e quations numeric ally using a 5th order 
Cash-Karp Runge-Kutta dPress et al.l [T992) method with adap- 
tive time stepping. As the electric and magnetic fields are pro- 
vided by the 3D numerical MHD experiment, these values are 
only known at the grid locations provided by the MHD experi- 
ment. Values at the actual particle positions are obtained using a 
Bi-cubic interpolation scheme. 

The data from the MHD simulations is saved in dimension- 
less units, and are scaled to SI units before being used in the trac- 
ing code. The electric field is not saved in the MHD code, and 
one can use different methods to derive the resistive contribution 
from the saved data Numerical codes typically uses resistivity 
models at constant 77, but with various forms of current depen- 
dent 77 threshold, higher than 2nd order current dependence - 
typically 4th order - to localise the effective rj towards the nu- 
merical resolution limit combined with a "shock" capturing ap- 
proach to localise diffusion in current sheet regions as long as 
new magnetic flux is advected into the sheets. In this case we 
have used the algorithm that was used in the numerical MHD 
experiment. The MHD code uses a high order finite difference 
approach on staggered grids to solve the MHD equations. The 
high order approach (6th order derivatives and 5th order interpo- 
lations) requires the special approach with 4th order dependence 
combined with the "shock" capturing mechanism to avoid/limit 
numerical ringing effects in the physical variables in the vicinity 
of regions where the magnetic field changes significantly over 
a few grid points. A descrip tion of the algorithm is given in 
iNordlu nd & Ga lsgaard! dl995h . 



The typical result is that the fastest particles reach energies much 
hig her than found in solar observations. 

lAmbrosiano et alj {1988) showed, using a 2D reconnection 
experiments, that having a single large scale diffusion region 
with internal structure gives rise to a stochastic acceleration pro- 
cess where initial close neighboring particles can experience 
very different accelerations histories as their paths diverge ex- 
ponentially with time. This is clearly different from earlier ex- 
periments that consisted of laminar diffusio n regions. The same 
process is seen in the previous 3D papers dTurkmani et al.l l2006: 
IDmitruk et al]|2003l 120041; IDmitruk & Matfhaeusll2006l) . where, 
as a general concept, the whole numerical domain is in a tur- 
bulent state with acceleration sites scattered randomly through- 
out the numerical domain. In the present case we have a situ- 
ation where the acceleration region is localised in the numeri- 
cal domain, and the global topology of the magnetic field is not 
in a turbulent state as in the earlier mentioned 3D experiments. 
This does not imply that the diffusion region has a simple lam- 
inar structure, as it contains significant variations in both space 
and time and we therefore expect the acceleration process to be 
stochastic in nature. Therefore, it is interesting to investigate how 
particles are accelerated in this environment, which type of par- 
ticle distribution arises, where they preferentially "exit" the free 
flight domain and which type of observational spectra are ob- 
tained from the bremsstrahlung generated as the particles impact 
with the photosphere. 

To initiate the investigation, a sub volume of the MHD box 
is chosen to contain the traced particles. The lower boundary 
of this domain is the photosphere in the MHD model, below 
which the plasma density becomes so high that collisions will 
thermalize any accelerated particles almost instantaneously. The 
top boundary is chosen such that for all the snapshots in the nu- 
merical experiment it is above the current sheet region that is 
responsible for accelerating the particles. The extent along and 
across the loop is chosen to cover the maximum extent of the 
current sheet. A representation of the particle box seen in rela- 
tion to the magnetic field configuration in the emergence phase 
is shown in Fig. [2] 

We perform tracing analysis for all snapshots provided by 
the numerical MHD emergence experiment. For each of the 160 
snapshots we initiate 90,000 particles in the chosen sub do- 
main. Inside this domain, they are given random positions and 
pitch angles. Finally they are assigned a velocity according to 
a Maxwellian distribution with a temperature of 1.2 MK. The 
particles are then traced for 0.5 seconds or until they either im- 
pact with the photosphere or exit the full MHD box above the 
photosphere. 

Tracing calculations are conducted for both protons and elec- 
trons. The results show that the parallel resistive electric field 
only provides for a small acceleration of the protons, while the 
much lighter electrons experience acceleration, even on multiple 
sites, that changes the distribution function significantly within 
the alloted tracing time. The following discussion will therefore 
only concern the results from the electron experiments. 



4. Particle setup 

Test particle tracing has previously been done in a number of 
different environments. Mostly thes e have consisted of situa- 
tions representing MHD turbulenc e ( Dmitr ulTet1ii|2003L r2004t 
IDmitruk & Matthaeus 2006; Turkman iet alJl2006h finding the 
accelerated particles to reach high veloc ities, having power law 
distributions with steep power indexes (Turkm ani et al.l l2006). 



5. Results 

In this section we discuss the results obtained from test particle 
tracing in the static MHD snapshots. There are three objectives 
to this section. The first is to follow the particle trajectories to see 
where they end up after the 0.5 seconds. As a fraction of these 
reach the photosphere, spatial distribution pattern indicates what 
we may expect to see in real observations, while the ones exiting 
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Fig. 2. The image represents the magnetic field line structure at 
one instance in the MHD simulations, showing the general struc- 
ture of the magnetic field by tracing a number of representative 
magnetic field lines. The dots represent the starting positions of 
a smaller sample of the particles in the tracing experiments. 

the top boundary are expected to enter into the solar wind parti- 
cle population. The second is to look at the velocity distribution 
functions obtained as a function of MHD snapshot to see how 
large a fraction of the particles are being accelerated and which 
energies they reach. Finally, we use the distribution function for 
the electrons that impact with the photosphere to derive spectral 
information that can be compared directly with observations. 

5.1. Impact patterns 

Following 90,000 particle paths through the MHD domain is by 
no means feasible, but provides detailed insight to the individual 
reasons for particle acceleration. We will look at a few examples 
of this below. Apart from this, there are two other approaches 
we are going to investigate in this section. One is to follow a 
sample of particles as they are moving around the domain. The 
most interesting are the ones that impact either the photospheric 
or the top boundary. These give us an idea about the locations 
where we would expect to see secondary effects of the recon- 
nection process, and allows us to obtain distribution functions 
and spectra of the impacting particles. 

5.1.1. Individual acceleration patterns 

In previous papers describing particle acceleration, it has 
been seen that the particles undergo repeated acceleration and 
de-acceleration as they pass through the turbulent magnetic 
field domain dDmitruk et al.i r2004; D mitruk & Matthaeusl l2006; 
iTurkmani et al.ll2006l) . In this experiment we assume that a sin- 
gle coherent current sheet is more likely to only provide one 
single acceleration of the particles, at least for the ones that very 
rapidly reach the domain boundaries. Fig.[3]shows the time evo- 
lution of a number of electrons. This shows a time history that 
is similar to those seen in the previous papers, starting with a 
rapid acceleration followed by a period of cyclic behavior in ki- 




Fig. 3. The graph represents the kinetic energy of a small num- 
ber of significantly accelerated particles as a function of time. 
The different colours represent different particles. The peak en- 
ergy in the diagram, 10 4 keV does not encompass the most en- 
ergetic particles in the experiment. It is noticed that all example 
particles are accelerated to energies many orders of magnitude 
above their thermal energies. About half of the sample particles 
stays at thermal energies. 

netic energy. We can conclude that there are a number of par- 
ticles that are trapped on trajectories where they repeatedly are 
moving forward and backward along the same magnetic field 
line. A small number of particles are seen to reach a peak energy 
level, which is subsequently represented by a line with constant 
energy. These are particles that have reached one of the imposed 
boundaries, from where further tracing is not possible. As can be 
seen, this type of particles reach the boundary rather fast, hav- 
ing experienced a rapid acceleration over a small fraction of a 
second. 

From inspection of the animations of the particle motions 
(movies are available with the Web version of the paper) it is 
easy to see both the fast moving particles that exit the domain 
and the ones that seem to oscillate around an equilibrium posi- 
tion. 

5.1 .2. The global view 

The behavior of the test particles depends strongly on the time of 
the snapshot from the numerical MHD experiment. Initially the 
particles are placed in a domain where no interaction between 
the two flux systems takes place. No parallel electric fields exist 
and the particles are just thermally moving along the magnetic 
field lines without significant changes in their energy. 

As the flux tube starts to rise into the region where the parti- 
cles are placed, effects on the particle behavior become visible. 
First the particles stay almost confined to the box in which they 
are initiated, but later they start moving around within this re- 
gion gaining an increasing amount of energy. As the tube clearly 
emerges above the photosphere and a strong current sheet forms, 
more action is seen by a fraction of the particles. In Fig. [4] two 
snapshots are shown. Here the box represents the full MHD do- 
main. The field lines indicate the structure of the magnetic field, 
while the colour of the particles indicates the energy obtained 
after 0.5 seconds of tracing. Looking at the images it is seen 
that only a small part of the volume becomes energized and that 
these particles are accelerated away from the reconnection re- 
gion following the global structure of the magnetic field. What 
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Fig. 4. The frames show the large MHD domain, with the im- 
posed particles at their final positions for two different times in 
the numerical MHD experiment. The field lines outline the gen- 
eral structure of the magnetic field. The particles are colored rel- 
ative to their final kinetic energy. 



is not clear from these images is the shape of the particle energy 
distribution function as they are channeled into the photosphere, 
where they are assumed to be de-accelerated almost instanta- 
neously, creating strong bremsstrahlung radiation (we will return 
to this below). 

It is clear from the images in Fig.|4]that only a small fraction 
of the particles are exposed to a very strong acceleration (larger 
than 1 MeV). To make this clear, one can derive an energy dis- 
tribution function after the 0.5 second integration. Doing this we 
ignore the fact that a large fraction of the particles that reach the 
defined boundaries do it well before the the 0.5 second has been 
reached. These are included as having the energy at the time of 
exit of the domain. Fig.[5]shows the time evolution of the kinetic 
energy distribution function. Initially we see the Maxwell distri- 
bution. As time in the numerical MHD experiment progresses a 
clear high-energy tail develops. Only minor adjustments occur 
between the different MHD snapshots. The power slope of the 
distribution function is roughly maintained, while the peak en- 
ergy changes, dependent on the level of activity in the snapshot. 
Looking at a single distribution function one can see different 
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Fig. 5. The figure shows the time evolution of the particle kinetic 
energy distribution as a function of time in the numerical MHD 
experiment. 



power slopes along the distribution, and it is possible to identify 
three different power indexes. 

The peak energy reached in the experiment is orders of mag- 
nitude larger than observed in sola r flares. This same o bservation 
is noticed in a previous analysis (Turkmani et al. 2006). We will 
return to this in the discussion below. 

5.1.3. Impact structures 

Direct observations of particle paths in the solar atmosphere is 
not possible. What is observed is the effect as the particles im- 
pact with the dense plasma in the photosphere. The radiation 
indicates where the field lines on which the accelerated parti- 
cles move intersect with the photosphere. From the tracing ex- 
periments it is simple to make images showing where and with 
which energy the particles exit the allowed tracing domain. The 
most direct comparison can be made in the photosphere, were 
fast particles will loose their kinetic energy in exchange for heat- 
ing the local plasma in a fraction of the observation time. This 
gives rise to local hot spots seen as H a ribbons in flares. We 
plot the positions and energies of the particles that reach the 
photosphere in Fig. [6] The frames show three snapshots repre- 
senting different times in the numerical MHD experiment. The 
colour represents the energies of the particles impacting with the 
photosphere at any time in the 0.5 second integration time. The 
numbers in the lower left of the frames indicate the number of 
particles (out of the total number) that impact with the photo- 
sphere. The frames show a number of interesting features. In the 
early phases only few and low energy particles impact. As the 
reconnection sets in a characteristic pattern develops, an almost 
closed envelope with a thick finger running down through part 
of the region. What does this signature express? And why do the 
later structures show much more complexity at especially one of 
the outer edges? These questions need to be answered in order 
to better understand the acceleration process. 

From t he topology analysis of the numeric al MHD experi- 
ment (See Galsgaard & Mor eno-Insertisl I2009L for details) it is 
found that the reconnection process creates three independent 
flux domains. One represents the initially open magnetic field re- 
gion, one confines the emerging magnetic flux and a final region 
represents the new loop system created by the ongoing magnetic 
reconnection process. Fig. Q] above shows images of how the 
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Fig. 6. The pattern and energy distribution of the particles after 
0.5 seconds time integration. The colour indicates the particle 
kinetic energy on impact with the photosphere. The three images 
represent different times in the numerical MHD experiment. 



the separatrix is a mathematical surface defining the boundary 
between two flux regimes. The intersection of the separatrix sur- 
face with the photosphere therefore defines a mathematical line. 
On the other hand, the current sheet, and the associated diffusion 
region, has a finite spatial extension. The field lines penetrat- 
ing the diffusion region therefore cover a much larger area when 
traced down to the photosphere. The spatial extent of the acceler- 
ated particles, when impacting with the photosphere, must have 
a finite extent centered on the separatrix surfaces intersection 
with the photosphere. 

The structure and location of the particle impact with the 
photosphere is therefore a direct consequence of the extent and 
mapping for field lines penetrating the diffusion region and con- 
necting to the photosphere. But, why does it include these small 
structural variations seen in the last of the images? Looking at 
the same type of impact map further up in the atmosphere, the 
mapping contains a much simpler mapping structure. The rea- 
son for the complicated patterns is located in the lower parts of 
the atmosphere. Here the magnetic fields lines are found to bend 
in different directions due to a weak convection like motion tak- 
ing place. This indicates that the pressure gradient locally domi- 
nates the momentum equation, implying that the plasma density 
may be high enough to thermalize the accelerated particles at 
this height in the atmosphere. 

From comparing the images in the two figures it is clear that 
the impact region maps directly back to the acceleration region. 
Being able to make realistic field extrapolations from magne- 
tograms into the corona, will therefore give a clear indication of 
where the diffusion region will be located in 3D space. Doing 
this with large precision is not a simple task in a magnetic field 
that is far from being in a potential or even linear force free state. 

5.2. Radiation spectra from the photospheric impact 

From a physical point of view the distribution functions give di- 
rect insight to the behavior of the particles. But, it is not a quan- 
tity that can be easily observed. Instead, we observe spectra from 
the radiation particles emit as they are de-accelerated on im- 
pact with the photosphere. This process occurs almost instanta- 
neously which allows us to ca lculate the spe ctra directly from the 
particle distribution function. Brownl dl97ll) and lHolmanl (2003) 
provide detailed information about how to do this. Because of 
the very high energy tail of the distribution function, the rela- 
tivistic formulae for the derivation must be used. This limits the 
energy of the most energetic photons in the spectrum, as it would 
otherwise only show one power law extending to very high en- 
ergies. 

The spectra obtained from this exercise is shown in Fig. [7] 
This shows two situations, both containing three different power 
law spectra depending on the energy range of the spectrum. The 
slope of the spectrum, covering the X-ray to gamma-ray range, 
changes from about 1.15 and 1.90. These slopes represent much 
harder spectra than observed from flare events. For large flares, 
the peak energy typically reaches hundreds of MeV and the X - 
ray power index is typically in the range of2.6to3.9 dLinl2008h . 



connectivity pattern of the numerical MHD experiment changes 
with time in the photosphere, as already discussed. 

Comparing Fig. Q]and Fig. [6] reveals close similarities be- 
tween the patterns of the two analyses. This is expected, as the 
region where particles are accelerated coincides with the loca- 
tion of the intersection of the separatrix surfaces. On one hand 



6. Discussion 

The previous section discusses the results we have obtained trac- 
ing test particles using MHD snapshots from a numerical flux 
emergence experiment. We have found that even in a situation 
where there is only one large current sheet located inside the nu- 
merical domain, is it possible to accelerate particles away from 
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Fig. 7. The spectra obtained directly from the electron kinetic 
energy distribution functions for the electrons impacting with the 
photosphere in the tracing experiments. The assumption is that 
the particles are thermalised over a short timescale (compared 
with the observation time) and distance through bremsstrahlung. 



the initially given Maxwell distribution. The final result is an 
electron distribution function that contains a high energy power 
law tail which extends from the Maxwell distribution. After the 
onset of MHD reconnection, this distribution forms rapidly and 
maintains a similar structure throughout, although the peak en- 
ergy and fraction of accelerated particles varies with time. There 
are a number of critical issues related to the approach used here: 

- Individual particles follow paths that gyrate around magnetic 
field lines. As part of the analysis we have compared the re- 
sults obtained using both the correct Lorentz equation for 
the particle motion and the much faster GCA. It is clear that 
for this experiment, the two methods produce almost identi- 
cal results. The reason for this comes from the fact that the 
background field from the numerical MHD experiment does 
not have a grid resolution that is even close to the Larmor ra- 
dius of the particle orbits. The magnetic and electric field 
therefore have spatial variations on length scales that are 
many orders of magnitude larger, which makes the GCA a 
good description of the particle motion. There is therefore 
no reason to use the LE description as long as the numeri- 
cal MHD experiment describes phenomena on typical solar 
surface length scales. 

- A limitation using the test particle tracing is the lack of in- 
teraction between the particles present in the domain. Each 
particle is treated independently of all other particles in the 
domain. This has a number of problems: 

- There are no built in collisions, which allows particles to 
be accelerated far more than in a realistic plasma. 

- There is no electric interaction, implying that the charge 
separation between the electrons and the much slower 
protons is again unrealistic. 

- Even if we included simple particle collisions, then a num- 
ber of particles would reach high energies and their collision 
cross sections would decrease sufficiently for them to escape 
the acceleration region. But, as the particles are accelerated, 
there is a natural exchange of energy with the macroscopic 
field. This would decrease the E-field in the diffusion region 
and the remaining particles would feel a smaller acceleration 
force, making it less likely that they would experience the 
same significant acceleration and escape from the region. 



- The above items have significant implications on the ob- 
tained distribution functions. In the present approach too 
many particles reach far too high energies. The lack of 
back reaction on the macroscopic fields allows the produc- 
tion of power laws that are too hard compared with ob- 
servations. These problems naturally influence the derived 
bremsstrahlung spectrum generated by the particles impact- 
ing with the photosphere. 

To improve on these issues we need a different and much 
more calculation heavy approach. One way is to adopt a PIC ap- 
proach, where one assumes the included particles to have a very 
exaggerated weight to allow the domain to cover the MHD scales 
without being dependent on resolving spatially the Larmor radii 
of real charged particles. This approach is being investigated. 

7. Conclusions 

The usage of test particle tracking in MHD simulations consti- 
tutes a simple way to obtain information about regions of particle 
acceleration and the spectra and spatial distribution on planes of 
interest. When the MHD data represents a large solar structure 
and the typical length scale involved is much larger than the par- 
ticle Larmor radius, the Guiding Center Approximation is the 
favored method to obtain these results. Comparisons between 
the GCA and the full Lorentz solution for this data discussed in 
this paper show no differences in results, while on the computa- 
tional side the GCA is between a 100 and a 1000 time cheaper 
in computing time. 

Warnings must be given. The lack of back reactions between 
the accelerated particles and the provided MHD field results in 
far too many particles being accelerated to far too high energies. 
The distribution functions and spectra obtained from such cal- 
culations are therefore giving too many high energy particles re- 
sulting in too hard spectra when calculating the bremsstrahlung 
emitted by the particles interacting with the dense photosphere. 

On the positive side, it is possible to use the impact loca- 
tion on the photosphere to identify the location of the separator 
surfaces connecting up to the diffusion region where the accel- 
eration is taking place. This information may be a help to guide 
magnetic field extrapolations of the photospheric field distribu- 
tion, as one can both see the impact points in the photosphere 
and the X-ray region in the corona where the reconnection pro- 
cess takes place. These have to fit with the topological structure 
obtained from the extrapolation models. 

To obtain more generally useful information from this type 
of exercise one has to adopt a more physically consistent model, 
which takes the whole approach into a much more complicated 
problem. But one we have to find a solvable solution to if we 
want to obtain reliable information from numerical experiments. 
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